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Abstract 

We study the stability of the BTZ black hole in the new massive gravity. This is a nontrivial 
task because the linearized equation around the BTZ black hole background is a fourth order 
differential equation. Away from the critical point of m^£^ = 1/2, this fourth order equation 
is split into two second order equations: one describes a massless graviton and the other is 
designed for a massive graviton, which could be obtained from the Fierz-Pauli action. In this 
case, calculating quasinormal modes leads to confirm the stability of the BTZ black hole. 
At the critical point, we derive two left and right logarithmic quasinormal modes from the 
logarithmic conformal field theory. Finally, we identify two s-massive modes propagating on 
the black hole background through the conventional black hole stability analysis. 
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1 Introduction 



It is well known that Einstein gravity in three dimensions has no propagating degrees of 
freedom. Massive generalizations of three-dimensional gravity allow propagating degrees 
of freedom. Topologically massive gravity (TMG) is the well-known example obtained by 
including a gravitational Chern-Simons term with coupling /i [H [2] . The model was extended 
by the addition of a cosmological constant term A = — to cosmological topologically 
massive gravity (CTMG) [3]. The gravitational Chern-Simons term is odd under parity 
and as a result, the theory shows a single massive propagating degree of freedom of a given 
helicity, whereas the other helicity mode remains massless. The single massive field is realized 
as a massive scalar = z^^'^hzz when using the Poincare coordinates and z covering the 
AdSs spacetimes [1]. However, it was claimed that the negative-energy massive graviton 
disappears at the critical point of /i£ = 1 [5]. This cosmological topological massive gravity 
at the critical point (CCTMG) may be described by the logarithmic conformal field theory 
(LCFT) [SI E] even for the zero central charge cl = 0. 

Another massive generalization of Einstein gravity in three dimensions was proposed re- 
cently by adding a specific quadratic curvature term to the Einstein-Hilbert action [HI |9]. 
This term was designed to reproduce the ghost-free Fierz-Pauli action for a massive propa- 
gating graviton in the linearized approximation. This gravity theory proposed by Bergshoeff, 
Hohm, and Townsend (BHT) became known as new massive gravity (NMG). Unlike the 
TMG, the NMG preserves parity. As a result, the gravitons acquire the same mass for both 
helicity states, indicating two massive propagating degrees of freedom. It was shown that 
there is no ghost in the linearized BHT gravity by performing a canonical analysis in fiat 
spacetimes [10]. Furthermore, in fiat and de Sitter spacetimes, the authors jllj have found 
two massive propagating scalars a and derived from the metric perturbations, satisfying the 
Klein-Gordon equation [V^ — m^]{cr, 0} = 0. However, up to now, there is no explicit form 
of two massive scalars propagating on AdSs spacetimes, even for the fourth order linearized 
equation for graviton was known under the transverse and traceless gauge |12] . 

It is well known that the BTZ black hole as solution to Einstein gravity with A is also 
a black hole solution to the CTMG. However, this does not necessarily imply that there is 
no difference in the dynamics of perturbations. It is obvious that perturbation discriminates 
between Einstein gravity and CTMG. Recently, it was shown that the (non-rotating) BTZ 
black hole is stable for all values of fi against the metric perturbations in the TMG by 
considering left-and right-moving normal modes ^13j. They have confirmed the stability by 
solving the massive scalar equation of (Vl^^ — m?)ip = 0. 

In this work, we wish to perform stability analysis on the BTZ black hole in the NMG, 
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which is a nontrivial task. The basic idea of performing the black hole stability is that 
one decouples the second order linearized equations and then, manages to arrive at the 
Schrodinger-type equation with an effective potential for physically propagating fields [T^ITS]. 
If all potentials are positive for whole range outside a event horizon, the black hole under the 
consideration is stable. If an effective potential is not positive definite everywhere outside 
a horizon, a special trick called the S'-deformation technique may be used to prove the 
stability [16] . It is well known that a practical tool for testing stability of all kinds of black 
holes is a numerical investigation of quasinormal frequencies u = ur — iui by imposing the 
boundary condition: ingoing waves near a event horizon and the Dirichlet boundary condition 
at infinity [T7]. That is, the unstable mode is defined by the condition of 

ijR = 0, ui<0. (1) 

in quasinormal mode approach [18]. We wish to perform the stability analysis of the BTZ 
black holes by computing quasinormal modes of the NMG. However, it seems that this 
analysis is not a straightforward task in the NMG because the linearized equation is a fourth 
order differential equation. If the mass parameter is off the critical value (m^ ^ l/2£^), 
the fourth order equation split into two second order equations with transverse and traceless 
gauge conditions: one is for a massless graviton (gauge artefact) and the other is for a 
massive graviton, which takes a similar form obtained from the Fierz-Pauli action. However, 
at the critical point of = l/2£^ (c^ = 0,c/j = 0), the fourth order equation leads to 
[V^ — 2A]'^h^i, = 0, which is difficult to be solved unless the LCFT is introduced. Hence, 
this case is surely beyond the standard stability analysis of a black hole prescribed above. 
Recently, Sachs and Solodukhin have determined quasinormal mode of black hole spectrum 
for tensor perturbations in the TMG [T9]. In their operator calculation, they have used 
the chiral highest weight condition of Lih^^ = to derive quasinormal frequencies, which 
are similar to the scalar quasinormal modes. However, this method is inappropriate to 
derive quasinormal modes at the critical point /i = l/£ of the TMG. To this end, Sachs has 
generalized the one-to one correspondence between quasinormal modes in the BTZ black hole 
and the poles of the retarded correlators in the boundary conformal field theory to include 
logarithmic operators [20] • On the other hand, Liu and Wang have studied the stability of 
the BTZ black string against gravitational perturbations in four dimensions, which is very 
similar to the Fierz-Pauli action in three dimensions |21] . 

The organization of our work is as follows. In section 2, we briefly review how the BTZ 
black hole comes out from the NMG. We derive quasinormal modes of the NMG by solving the 
first order equations with A = — 1, which is based on Sachs and Solodukhin method used in 
the TMG in section 3. However, the left- and right-moving modes are not orthogonal, which 
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has a problem to be considered as two independent massive modes. Section 4 is focused on 
deriving quasinormal modes at the critical point of = 1/2 with A = — 1. We identify two 
massive modes $ and from the NMG (Fierz-Pauli action) using the conventional black 
hole stability analysis in section 5. Finally, we discuss similarity and difference in stability 
between the TMG and the NMG in section 6. 

We would like to mention the choice of a cosmological constant A. In general, we choose 
A = except section 3 and 4, where it is chosen to be A = —1 for the simplicity of 

operator computation. 



2 New massive gravity 

The NMG action [S] composed of the Einstein-Hilbert action with a cosmological constant A 
and higher order curvature terms is given by 

q(3) _ q(3) I q(3) /r)\ 

'-'NMG ~ '-'EH "T ^//C K'^J 
1 

IQttG 



S^H = I d\^g{R-2\l (3) 



4^ - -T^/^^-v/^(^-^'^^-^^> 

where G is a three-dimensional Newton's constant and m? a parameter with mass dimension 
2. From now on, we set G = 1/8 for simplicity. The Einstein equation is given by 

1 

where the Einstein tensor G^j, and K^j^y tensor are given by 



3 3 
+ 4:RfipuaR^'^ — -RR^iu — RpuR''" g^u + -R g^u- (6) 

Z o 



In order to have a black hole solution with dynamical exponent z [221 123] , it is convenient to 
introduce dimensionless parameters 



y = m' w = X r, (7) 



where y and w are proposed to take 



2^-3^ + 1 + Z + 1 

y = t; > ^ = ^ — • 
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For z = 1 (nonrotating) BTZ black hole, one has y = ^ and w = — |, while y = —\ and 
w = — Y are chosen for z = 3 Lifshitz black hole. 

In this work, we consider the BTZ black hole solution only 



ds'^Tz = Qf.udx^dx'' = - \^-M + ^ J de + + (9) 

where is the ADM mass determined to be = with r+ the horizon radius. Impor- 
tantly, the mass parameter th? and cosmological parameter A are fixed as 

to obtain the BTZ black hole. In this background, taking into account g^'^K^j^ = 2m^A, the 
trace of ([5]) leads to the constant curvature scalar as 

i? = 4A = -|, (11) 

which is the same form as in the Einstein gravity {R = 6A) with the cosmological constant 
A = —l/i"^. On the other hand, the Ricci tensor takes the form 

1 4 _ 

Rfiu = \g^,u + ^^^f^'^ = ^^9fiu, (12) 
which is the same as that of the Einstein gravity 

R^, = 2Ag^,. (13) 
The curvature tensor R^pya takes the form 

R^,pya = A (g^^ugpa " gpaQpu^ ■ (14) 

For an AdS-sized black hole with r+ = i, one chooses the unit mass of = 1, which is 
designed for finding its quasinormal modes (frequences). Setting £^ = 1 (A = — 1), we rewrite 
the line element in global coordinates: 

(is^^i = — sinh^(p)(ir^ + cosh^(p)(i0^ + dp"^, (15) 

where the event horizon is located at p = (r_|_ = 1) while the infinity is at p = oo (r = oo). 
The black hole temperature is = 1/Att. The metric tensor g^^, can be when using the light 
cone coordinates u/v = t ± (f) 



( \ -icosh(2p) \ 

-icosh(2p) \ 

V i; 



(16) 
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Then the metric tensor (ITB]) admits the KiUing vector fields L^, k = 0, —1, 1 for local 
SL(2,i?)xSL(2,i?) algebra as 



cosh(2p) 



du 



1 



sinh(2p)^" sinh(2p)^''^ 2^^^. 



(17) 



and Lq, L_iji similarly by substituting u ^ v. Locally, they form a basis of the SL(2, R) Lie 
algebra as 

[Lo,L±i] = tL±i, = 2Lo, (18) 

which are useful for generating the whole tower of quasinormal modes. 



3 Quasinormal modes 

Considering the perturbation /i^j, around the BTZ black hole background in (fT6i) 

gfiu = g^iu + h^y, (19) 
the linearized equation to ([5]) takes the form 



5Gf,y{h) + \h 



1 



2m? 



5K,,{h) = 0, 



(20) 



where 5G^y{h) and 5K^y{h) are the linearized Einstein tensor and linearized K^y. We choose 
the transverse and traceless (TT) gauge to find a massive graviton propagation on the BTZ 
black hole background as 

V^/i'^^ = 0, h = hpP = Q. (21) 

We wish to mention that this covariant gauge is also convenient for studying the conventional 
stability of black holes in section 5. Then, the fourth order equation fl2U]) is split into 



2A 



2 5A\ 

H 

2 



h 



0, 



(22) 



which implies two branches of solutions [12]. In this section we choose A = — 1 for simplicity 
of operator calculation. The first equation is 



2A 



/if;/^ = 0, 



(23) 



whose solution corresponds to the unphysical modes of left- and right-moving massless gravi- 
tons, while the second is the equation describing a physically massive graviton with 2 DOF 



2 5A\ 

H 

2 \ 



h 



M 



0. 



(24) 



In order to have a non- negative mass, one requires > 1/2. For = 1/2 (the BTZ black 
hole background), the fourth order equation (1221) could not be split into the two second order 
equations ([23]) and (EH). 

First of all, we consider the AdSs background [or M = — 1 in ([9])] 

'^■^AdSa ~ ~" cosh^(p)(ir^ + sinh^(p)(i0^ + dp^^ (25) 

one can easily find the solution [5] to (124|) as 



vM,AdS3 
"'Hu 



-ihu—ihv 



sinh^(p) 
(cosh(p))^+^ 





1 


h-h 
2 


ia 


\ 




h-h 
2 


1 


ib 




v 


ia 


ib 


-a" 


1 



(26) 



where a and b are given by 



h-h 



sinh(2p) ' sinh(2p) 
We have two sets for h and h as primary states 

2 + V2 + 4m2 



h + h 



h-h = ±2, 



(27) 



(2J 



where ± denote the solution to the first order equations {D'^^ ^ '^^ h) = with the operators 
in Eq. f l30|) . respectively. This implies that the solutions to the first order equations are also 
those to the second order equation. In addition, the solution to fl23|) in the AdSs background 
(I2S]) has the same form as in ( 12^ when substituting {h, h) = (2, 0) for left-moving mode and 
(0, 2) for right-moving mode: these are also the solutions to {D^/^h^/^)^^ = 0, respectively. 

In this sense, hereafter, we use mainly the first order equations instead of higher order 
equations to find quasinormal modes. The fourth order equation (122|) can be expressed 



in terms of mutually commuting operators as 







with 



M 



m? + 



(29) 



(30) 



(31) 



At the critical point (the BTZ black hole background: m? = 1/2, M = 1), the operators -D*^ 
and degenerate, while the operators D^^ and degenerate. Thus, we have to treat it 
separately in the next section. 
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On the other hand, the second order massive equation can be expressed as 



(32) 



In order to find quasinormal modes, we have to solve (l32l) together with the TT gauge (l2Ti) . 
This approach was used to derive quasinormal modes in the TMG [19], which shows the 
quasinormal modes of a minimally coupled scalar. We remind the reader that the solution 
to fl32p may be a linear combination of two solutions to the first order equations 



O^e. "'^Vahsu + Mh 



/ii/ 



^ e, 



a/3 



Mh 



0, 




(33) 
(34) 



together with the TT gauge fl2Tl) . Let us first find the right-moving solution to fl33|) . The 
least damped {n = 0) quasinormal mode to fl33l) could be found by considering the form 



h 



M 



where 







i 

sinh(2p) 





\ 



(35) 



(36) 



\ ^ 3 / 

\ sinh(2p) " sinh2(2p) / 

The transversality condition of V^/i^^ = implies the chiral highest weight condition of 
Lih^i, = under the form of hff^^ in (l35i) . giving the constraint 



2ip+ + 2ip_ cosh(2p) + sinh(2p)9p F*^(p) = 0, 
whose solution is given by 

F^^(p) = [sinh(p)]-2^P-[tanh(p)]-'^ 
From (l33l) for fi = u = p, p_ is determined to be 

p_ = -ihniM), hR{M) -- 
At the first sight, the n = quasinormal mode seems to be 



M 



1 

2' 



-ik{T+4,)-2hR{M)T^M (^PY 



Considering the form of quasinormal frequency 



(37) 



(3J 



(39) 



(40) 



(41) 
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we read off it from ( I4UI) 



ojm = k — 2ihji{M). 



(42) 



We note here that "fc" is the angular quantum number as well as the real part of quasinormal 
frequency. In order that um be a quasinormal mode, hn^M) is required to be positive and 
thus, one has to have M < — 1. However, it is impossible to make M < — 1 because M > 1, as 
is shown in f l3ip . Hence it is clear that fl4Up is not considered as a right-moving quasinormal 
mode for the least damped case. 

The next step is to find the right-moving solution to fl34l) . This solution is obtained simply 
by replacing M hj M = -M [equivalently, /ii?(M) by hn{M)] in as 



(43) 



where 



hn{M) 



M > 1. 



(44) 



M 1 
Y~ 2' 

It seems that fH5]) is regarded as a right-moving quasinormal mode of the least damped case 
whose quasinormal frequency is given by 



M 



k - 2ihR{M). 



(45) 



Hence, we could construct the overtone quasinormal modes for this solution. Introducing 
the operator combination which replaces uj hj uj — 2, the n^'^-overtone quasinormal 

mode is constructed by 

hf/' = {L-^L-^r hf^ (46) 
whose quasinormal frequency takes the form 



CO 



M 



k-2i 



hii{M) +n 



neZ. 



(47) 



Similarly, when imposing the anti-chiral highest weight condition of Likf^^ = and thus 



= —ih^lM), the left-moving quasinormal modes to (!33|) take the form 
where 



(4J 



^S(p) = [sinh(p)]-''^^^'"^[tanh(p)] 



ik 



/o 







1 

2 



\ 



sinh(2p) 

4 



hUM) 



M 1 



(49) 



sinh(2p) sinh2(2p) 
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whose quasinormal modes are given by 



M 



-k - 2i 



hL{M)+n 



ne Z. 



(50) 



We note that two quasinormal modes h'^l and /i^ are ingoing (the right- and left-moving) 



modes at the horizon p = and are normalizable modes at p = oo satisfying the Dirichlet 
boundary condition. It is worth noting that the left-moving quasinormal modes to with 



M 

T 



(51) 



is not available since hL{M) < due to M > 1. We summarize the possible quasinormal 
modes in the following Table. In the BTZ black hole background, (133|) allows the L-mode 



equation right-moving 


I (R) mode 


left-moving 


(L) mode 


m Ki'^ 


(X) 
(0) 


"'flV 


(0) 
(X) 



only, while f lM|) does the R-mode. 

In this section, we have obtained the L/R- moving quasinormal modes, which show at least 
the s-mode stability using ([1]): for A; = (wr = 0), one has a;"^^^^^^ = 2n + hji/L^M/M) > 0. 
In this case, k = requires the disappearance of the L/R-moving modes. 

At this stage, it seems appropriate to comment on the similarity and difference between 
TMG and NMG. In deriving the quasinormal modes, we have used the operator method. As 
far as this method is concerned, there is no essential difference between TMG and NMG. The 
difference is that the substitution of M — )■ p in fl33p is necessary for the TMG and f l34p is 
required additionally for the NMG because the parity is preserved in the NMG. However, we 
note that the L/R-moving modes f H6|) and fjlHj) are not orthogonal to each other even s-mode 
{k = 0) in the NMG because they contain hpp commonly. This gives rise to a significant 
difference in obtaining the quasinormal modes of the same BTZ black hole between two 
massive gravity theories. 

Finally, at the chiral point of M = 1, two /i^ and h^l are no longer n = quasinormal 
modes because of = (/^/^/^(l) = 0). Hence, we will treat it in the next section. In 
addition, we mention that requiring both of L^ih^^ = and L^ih^^ = lead to normal 
modes with ingoing and outgoing fluxes at the event horizon as in Ref. [13], which is not the 
condition for obtaining quasinormal modes. 
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4 Quasinormal modes at the critical point 



In order to see what happens in the perturbation at the critical point (the BTZ black 
hole background), we need to solve the fourth order linearized equation. In this section 
we choose A = —1 for simplicity of operator computation. In general, a mode annihilated 
by D'^ (D^)[D^]{(D^y but not by D^} is called massive (left-moving) [right-moving] {left- 
logarithmic (L,new)} and is denoted by hffj^{h^j^)[h^^,]{h^i!^^'^} [25]. Also, a mode annihilated 
by [D^Y but not by is called right-logarithmic mode and is denoted by h^i^*^^ . Away from 
the critical point {m? ^ 1/2), the general solution to fl29|) is obtained by linearly combining 
left, right, and two massive modes. 

At the critical point, rn? = 1/2, degenerates into D^, while D^^ degenerates into D^. 
The L/R- moving modes are purely gauge degrees of freedom (unphysical modes), whereas 
two massive modes and their logarithmic partners constitute physically propagating bulk 
modes. At the critical point (the BTZ black hole background), the fourth order equation 
becomes 



■\L\2 j^L / R,new 



^ 



2A 



j^L / R,new 



0, 



which is basically different from the second order equations fl23l) and (E 
when observing 



(52) 

The reason is clear 



//ii/ 



2A 



j^L / R,new 



(53) 



which was derived in Appendix I in detail. See Ref. [26j for the derivation this relation on 
AdSs background. This naturally leads to flS^ when operating (V^ — 2A) on both sides and 
using (123|) . In this case, considering M — )■ L and M — )■ i?, one has 



h 



Jk{T-<l>) 



[tanh(p)] 



ik 



[0 







1 

2 





2 



sinh(2p) 
4 



(54) 



sinh(2p) sinh^(2p) 



''fj.u 



~ik{T+4,)^lR 



-ik{T+4i) 



[tanh(p)] 



-ik 



( 



\ 





2 



sinh(2p) 









sinh(2p) 



4 

smh^(2p) 



(55) 



Then, we construct two newly logarithmic solutions of L/R modes to 



h^R,new 

"'fJ.L' 



M I 



dm.2 h 



a2=l/2 
1^=1/2 



fii/i 



(56) 
(57) 
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where yij^p) is defined by 



Z/(^> p) 



ln[sinh(p)]. 



(5J 



We are now in a position to construct two quasinormal modes in the BTZ black hole 
background inspired by left-logarithmic sector of TMG. Unfortunately, h^'J^'^^ is growing in 
time r and p, showing disqualification as a quasinormal mode. Hence, Sachs [20] has proposed 
that the left-logarithmic quasinormal mode can be constructed by considering 



"'fj.u 



vL,new 



For example, a candidate for quasinormal mode is the first descendent given by 



(59) 



+ 



^ -y(r, p) + yy(r, p) 



2e 



-2t 



sinh (p) 



cosh^(p) + -k'^yir, p)+ik[ y{T, p) 



1 + cosh^(p) 



(60) 



as the left-logarithmic solution. Here 



2e 



-2t 



sinh (p) 



(tanh(p)) 



ik 





1 



Y smh(2p) 



2 cosh(p) 
sinh(p) 



sinh(2p) 
2 cosh(p) 
sinh(p) 
4(1+2 cosh(2p)) 
sinh^{2p) 



(61) 



where we observe that hjl'}^^ shows a genuine quasinormal mode with exponential fall-off in 
r and p. Note that in the limit of A; — )■ 0, the first descendent becomes (2.18) of Ref. [20] for 
the TMG without ipf;^ which may be considered as as pure gauge □. Here we emphasize again 
that the A; — > limit denotes both the s-mode of (p and ur = 0, which confirms the s-left 
logarithmic mode stability of BTZ black hole. The higher descendent modes in principle 
can be obtained by using fl59|) . however, it seems to be a formidable task. Thus, we confine 
ourselves to the k = (s-mode) case, for simplicity, and obtain the second descendent mode 
as 



"'fii/ 



h 



L,new 



-At 



sinh p 



L,{2A) 



12?/(r, p)il) 



L,{2B) 



(62) 



*We thank I. Sachs for pointing out it. 
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where 



2(11+15 cosh(2p)) 
sinh(2/3) 

4(7 + 3cosh(2p))coth(p) 



2(11+15 cosh(2p)) 4(7 + 3eosh(2p))coth(p) M2e+'^^ cos1^p)+9 cosH^p)) 



( 6 ll + 9cosh(2p) 

11 + 9 cosh(2p) 11 + 9 cosh(2p) 

f 15 cosh 
\ sinh(2p) 

/ 1 2 + cosh(2p) 4coth(p) 

2 + cosh(2p) 2 + cosh(2p) (5 + cosh(2p)) coth(p) 

y 4coth(p) (5 + cosh(2p)) coth(p) 



\ 



2(3+2 cosh(2p)) 



(63) 



sinh'^ (p) 



Restoring the /c-dependence, we may find by induction that the higher modes h 



R,(n) 



'n>2) 

are constructed by replacing e"^"^ by e"^""^. Then, the quasinormal frequencies may be given 
by 



—k — 2m, n & Z, 



(64) 



which are the same quasinormal frequencies of left-moving modes with hR{M 



in dSD. 

The right- logarithmic quasinormal modes could be similarly constructed by 



hR,W 



pf 



+ 



2e 



-2t 



sinh (p) 



cosh^(p) + -k'^y{r,p) - ik ( ?/(r,p) 



1 + cosh^(p) 



h 



R 



where 



2e 



-2t 



sinh (p) 



-^'=("+'^)(tanh(p))- 



ik 



1 





2 cosh(p) \ 
sinh(p) 
2 



sinh(2p) 

2cosh(p) 2 4(1+2 cosh(2p)) 

y sinh(p) sinh(2p) sinh^(2p) / 



(65) 



(66) 



R (1) 

Here we observe that h^!/ is a genuine quasinormal mode with exponential fall-off in r and 
p. Considering the k = (s-mode) case, for simplicity, the second descendent mode is given 
by 



"-pi/ 



2 p- 
j^R,new 



-At 



sinh p 



2^5:(^^)-12y(r,p)^f'(^^) 



pi/ 



(67) 
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where 



R,{2A) 



( 10 + 7cosh(2p) 
9 + 3 cosh(p) + 5 cosh(2p) 



31+33 cosh(2p)+4 cosh(4p) 
sinh(2p) 



9 + 3 COSh(p) + 5 COsh(2p) 31+33cosh(2p)+4cosh(4p) ^^ 

21+25 cosh(2p) 
sinh(2p) 

21+25 cosh(2p) 2P(p) 
sinh{2p) sinh^(2p) 
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2 + cosh(2p) 
2 + cosh(2p) 
y (5 + cosh(2p)) coth(p) 



2 + cosh(2p) (5 + cosh(2p)) coth(p) \ 
1 4 coth(p) 

4 coth(p) 



2(3+2 cosh(2p)) 



(68) 



sinh'^ (p) 



with 



P(p) = 36 + 4cosh(p) + 40cosh(2p) + 10cosh(4p) + 6sinh(2p)) + 3sinh(4p). (69) 

Assuming the fc-dependence again, the higher modes h^u^^ {n > 2) can be obtained by 
replacing e"^"^ by e~'^^'^ as welL Thus, the quasinormal frequencies may be given by 



cu^ = k — 2m, n E Z, 



(70) 







which are the same quasinormal frequencies of right-moving modes uj^j with hL{M = 1] 
in ([50]). 

In this section, we have obtained L/R-logarithmic quasinormal modes, which show at 
least the s-mode stability using ([T]): for /c = {cor = 0), one has oj^ j^/ji = 2n > 0. 

Finally, we note the similarity and difference between TMG and NMG. In deriving the 
quasinormal modes at the critical point, we have used the operator method. There is no 
essential difference between TMG and NMG. The difference is that the L-logarithmic mode 
approach is necessary for the TMG, while both of L/R logarithmic approaches are required 
for the NMG because the parity is preserved in the NMG. However, it is not clear that 
two L/R-logarithmic quasinormal modes are orthogonal to each other even for the s-mode 
(A; = 0). This is an important issue because the orthogonality may guarantee two independent 
massive quasinormal modes at the critical point. 



5 Two s-massive modes in NMG 

First of all, we point out a few problems on the previous tensor quasinormal modes. Their 
frequencies were constructed as those of a minimally coupled scalar mode. The L/R modes 
of fl46p and fl48p are not orthogonal to each other even for s-mode {k = 1) because they 
contain hpp commonly. Hence, we did not confirm that the L/R- moving modes are two truly 
propagating modes in the new massive gravity. This asks directly what is the explicit form 
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of two physically massive modes propagating in the BTZ black hole. In this section, we wish 
to identify two massive graviton modes propagating in the BTZ black hole background. 

In order to perform the stability analysis of massive graviton with 2 DOF, we first recall 
Eq. ([9]) as the BTZ black hole background. In this background the perturbation equations 
fl24|) can be rewritten as 



2 5A 
m +- 



0. 



(71) 



This equation is similar to the Fierz-Pauli massive equation 



h 



+ 2A 



(72) 



which is derived from Fierz-Pauli action together with the TT gauge. In order to solve (17T]) 
explicitly with the TT gauge, we consider the following two distinct perturbing metric ansatz: 
the type I has two off-diagonal components and hi 



/ii/ 



/ ho{r) \ 

hi{r) 
\ ho{r) h{r) / 



(73) 



while for the type II, the metric tensor takes the form with four components Hq, Hi, H2, 
and as 



H,{r) Hi{r) \ 
Hi{r) H2{r) 
H^{r) j 



(74) 



At this stage, we note that the above two choices of (1731) and (17^ are working only for 
s-mode {k = 0) perturbation. Also, two are orthogonal to each other. In Appendix II, we 
have shown that k = case leads to two consistent decoupling processes for I and II metric 
perturbations, when solving fITB . 

Substituting Eq. fl73|) into Eq. flTTj) and eliminating hi from {t, 0) and (r, (p) components 
of flTTj). we obtain 



{{m^ - l/2f){ryf -M) + col} h^ + 



r^/p + M 2 rye - M 



+ - 1/2^ + 



(m^ - l/2e) h[ 



ho = 0. 



(75) 
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Figure 1: V$ graphs as function of r with fixed values £^ = 1, Ai = 1, 0; = !, and m = 10. 
For ^ l/2£^ case, V$ blows up at r = rp = ^yPM + u'^/{m? — l/2^2)(~ 1.005, near the 
event horizon), while for the critical = l/2£^ case, there is no the blowing-up point. 

Note that we consider only s-mode {k = 0) for components (t, r),(r, r) of Eq. fl7T|) . At the 
horizon r = r+ and the spatial infinity r = oo, a solution to the above equation behaves as 



/lo ~ (r - 



ln[r— r+] 



and ho r 



(76) 



where r+ = yjAi^. Introducing the tortoise coordinate r* in [dr* = dr/{—Ai + r"^ and 
redefining ujh = —iu and a new field $ as 

ho 



with 



g{r) = {{m^ - l/2t'^)r^/P - (m^ - 1/2&)M + col} 
Eq.( 17S]) can be written as the Schrodinger-type equation 

^ + [a;2-\4(c.,r)]$ = 0, 
where the a;-dependent potential takes the form 



(77) 
(78) 

(79) 



2 13 3M 
+ — - — - + 



3(m2 - l/2£2)V2(rV£2 - M) 



2(m2- 1/2£2)(27W -3rV£2) 



0) 



P {(m2 - l/2£2)(->l + ryp) - 

It is important to note that for > 1/2^^^ the above potential is always positive for whole 
range of r+ < r* < cxd even though it blows up at r = (see Fig.l), which implies that the 
type I-perturbation is stable. 
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On the other hand, plugging Eq. (174|) into Eq.l FTTj) and after manipulations we obtain the 
second differential equation for Hi{r) as 

{(m^ - l/2^2)(rV£2 - M) + r^/t - 2M/e + ul) H'l + ( ^fJ^l'^.A 



,2 1 /o/?2 , 



In this case, we also consider only s-mode {k = 0) for component (t, 0) of Eq. (17T|) . We have 
wave forms near the horizon and spatial infinity 

/7i~(r-r+)"'^vi755=e'"'^vi75^""t'-''+l, and H, ~ ^-2±V^2-2+i/2_ (33) 
Introducing the tortoise coordinate r*, and redefining Uh = —ioj and a new field \l/ 

with 

, JMim^ - 1/2&) + 2M/(i^ - (m^ - l/2&)ryP - rVi^ - oj^ , , 

/(^) = ^ =: — (84) 
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Eq.( 15T]) can be written as the Schrodinger-type equation 



+ K-Kt(w,r)]^ = 0, 



^5) 



where the a;-dependent potential is given by 



rn? + 



37W 



+ 



3r2(m2 + l/2£2)2(r2/£2 _ 



Aim' + l/2e){r^/e - M) 



& {Mm? + ?>M /2£2 - m2r2 /P - r2 /2P + a;2) 



6) 



It is important to note that for > l/2£^, the above potential is always positive for whole 
range of r+ < r < oo (— oo < r* < 0) (see Fig.2) even though it blows up at r = rp, 
which implies that the type Il-perturbation is stable. This has been confirmed by computing 
quasinormal modes numerically when using ( IHTj) in Ref. [21], instead of (185|) . 

Finally, we have shown that the positivity of two potentials leads to the stability condition 
of the BTZ black hole 

(87) 



2P' 



which is consistent with the condition of M > 



6 Discussions 



We have performed the stability analysis of the BTZ black hole in the NMG. We have derived 
quasinormal modes of the NMG by solving the first order equations, which is based on Sachs 
and Solodukhin method used in the TMG. However, we have observed that the left- and 
right-modes are not orthogonal, which has a problem to be considered as two independent 
massive modes. Hence, we have identified two massive modes from the NMG (Fierz-Pauli 
action) using the conventional black hole stability analysis. Furthermore, we have obtained 
left- and right logarithmic quasinormal modes at the critical point of m^£^ = 1/2 using the 
LCFT. 

We discuss similarity and difference in stability between the TMG and the NMG. In 
the case of the TMG, the BTZ black hole is stable for any Chern-Simons coupling constant 
/i > by computing (l33l) with replacing M by /i, while the stability of the BTZ black hole 
is guaranteed for m^£^ > 1/2 (M > t) by computing two first order equations (!33|) and 
In the TMG approach, the authors [13] have confirmed the stability by using that a single 
massive scalar of = z^hzz satisfies the Klein-Gordon equation in the AdSa background. On 
the other hand, the two left- and right-modes of the NMG [ fj46|) and fHHl) ] are not orthogonal 
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to each other, implying that the two are not considered as two independent massive modes 
propagating on the BTZ black hole background. This requires a reanalysis of stability of the 
BTZ black hole in the conventional black hole approach. We have obtained two propagating 
s-massive modes $ and ^ by solving the Schrodinger equations fl79l) and f l85|) directly, which 
are surely not the Schrodinger equations for a massive minimally coupled scalar. 

However, the difference between the NMG and the Fierz-Pauli action is just the presence 
of the critical point of m^£^ = 1/2 (cl = 0,cr = 0), whose equation is given by the fourth 
order linearized equation (152|) . This point exactly coincides with providing the BTZ black 
hole solution. Actually, there is no definite way to solve this equation. In the strict sense of 
stability of the BTZ black hole in the NMG, we have confined ourselves to solving (1521) by 
using the LCFT technique in section 5. However, this method is very restrictive and thus, 
there is no definite way to confirm its results at present. Particularly, the s-mode {k = 0) 
computation could be easily done to derive the left- and right logarithmic quasinormal modes. 

At this stage, we would like to comment on how we can apply the stability condition (I57|) 
of the BTZ black hole in the NMG to the dual CFT whose central charge is given by [Bl Ell EZ] 



It is well known that the exact agreement is found between the quasinormal frequencies of 
the BTZ black hole and the location of the poles of the retarded correlation function of the 
corresponding perturbations in the CFT f2E\. This has provided a confirmed test of the 
AdS3/CFT2 correspondence. If this correspondence is still valid for the NMG, one observes 
from (l88l) that the zero central charge appears at the critical point, while the unitarity of the 
CFT2 would demand the bound of > l/2£^ due to c^/zj > 0. For the non-critical case, we 
find the stability condition of M > i from (!87|) . which turns out to be the unitarity condition 
for the CFT2. In this case, we expect that the quasinormal modes obtained here could be 
derived from the location of the poles of the retarded correlation function of the corresponding 
perturbations in the CFT2. On the other hand, it seems that the quasinomal modes obtained 
at the critical point of M = i could be found from the retarded correlation function of the 
corresponding perturbations in the LCFT as was shown in the TMG ^U\. Hence, it is clear 
that the stability condition of BTZ black hole provides the unitarity condition of the dual 



Finally, we would like to mention that the linearized higher dimensional critical gravities 
were recently investigated in the AdS spacetimes [29] but their quasinormal modes are not 
studied in the AdS-black hole background. The non-unitarity issue of log gravity is not still 
resolved, indicating that the log gravity suffers from the ghost problem. 

Consequently, we have performed the stability analysis of the BTZ black hole in the NMG. 



cl/r 




2m2£2j- 



1 1 



(88) 



(L)CFT2. 
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It seems that the BTZ black hole is stable against the metric perturbation by computing 
quasinormal modes and observing two potentials. However, the stability at the critical point 
is not still completely proved because two s-modes [f l60|l and fl65|) ] and [ fl62|) and f l67|) ] are 
not orthogonal to each other, respectively. 
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Appendix I: LCFT approach to the critical point of NMG 

In this appendix, we will prove the equations of motion (153|) at the critical point in the BTZ 
black hole of the NMG using the LCFT approach. We first note that 

Loy = Loy = ^, Liy = Liy = 0, (89) 



and 



Lo/iJ, = tkh^,, Loh^, = 0. (90) 
Acting on the R-logarithmic mode of h^^"^^ with Lq{Lq), we have 

which show that h^^^"^ is not an eigenstate of Lq or Lq as like the new mode at the chiral 
point in AdSs of CTMG |6j, but it could be an eigenstate of the subtraction operator Lq — Lq. 
The two representations of Lq and Lq are now given by the Jordan cell form 



L 







(i^R,new \ / n 1 \ / uR,new \ 

Thus, it indicates that h^^^^ is a L-logarithmic partner of /i^^,. This shows again h^i^^^ (the 
logarithmic partner of /i^^^) and h^^ are an eignestate of the subtraction operator Lq — Lq: 



Lq — Lq 



l^R,new \ / 1 n \ / UR,'new 



Now, using the SL(2,/2) quadratic Casimir of = i(LiL_i + L_iLi) — Lq, the equation 
(133]) is written by 

{D^D^h^'''^'")^^ = -(V^ + 2)/;.^'^^"' = [2(L2 + P) + 4] /ij;,"^"^. (94) 
Making use of the following equations 

Lih^;;'"" = -(2 - zfc)e" tanh(p)/i2:,"^"', Ii/iJ;,"""' = 0, (95) 

we arrive at 

(Z^^D^/i^'"'^"')^, = -(V' + 2)/i2:,"^"' = -2h^,. (96) 
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Operating (V^— 2A) on with A = — 1 and using (12^ on the R-moving mode /i^^,, this leads 
to dSg): ^ 0, (^D^D^)2hR'^ew ^ Q remind the reader that Lih^i''^'"{= 0) 

in fl95|) is the "chiral" highest weight condition, but not Lih^^^^ 7^ as shown in fl95|) . This 
differs clearly from the CTMG at the chiral point of the AdSa background [6J, where one 
should impose the "highest weight" conditions of both Lih^i^'^'^ = Lih^i^'^'^{= 0) [7]. 

On the other hand, for the new L- logarithmic mode (jSH]), we follow the same steps as did 
in the new R-moving mode. In addition to (l89l) . we have 

Lo/iJ, = 0, Lohfi^ = -tkh^:,. (97) 

Therefore, the new L-logarithmic mode of hj^H^^^ satisfies 



Similarly, the two representations of Lq and Lq take the compact matrix forms 



l^L,new \ 

"'^lU 













l^L,new \ 




^ -ik 






{ 






(100) 



l^L,new 

^° ( '\ } ^ [ )[^ 1 ■ 

This shows again h^i^'^^ (the logarithmic partner of h^^) and h^^^ are an eignestate of the 
subtraction operator Lq — Lq: 

Now, making use of the following equations 

Li/^Jr" = 0, Ii/^Jr" = -(2 + ^k)e'' tanh(p)i/(r, p)/ij„ (101) 

where Li/t^^"^'"(= 0) is the "anti-chiral" highest weight condition, but not for Li/i^^f^"^ 7^ 0, 
we finally arrive at 

(D^D^/i^'"^")^, = -(V^ + 2)/iJ;;"^" = -2/iJ„ (102) 
which clearly confirms flS^ . 
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Appendix II: Full perturbation analysis in BTZ black 
hole background 

Considering the full /i^,^ components given by 

/ Httir) HtJr) HtJr) \ 



Hrtir) Hrrir) Hr^ir) 
\ H^t{r) H^r{r) H^^{r) J 



the tensor perturbation equation of Vbtz^hi' ~ + '^^j^/i;/ = l^ad to 



(103) 



{t, t); -2r\M + Ar^K - - A^r^)K - {{2Mk'' + 2AmV2) - Kr\?,M - Ar") 

-^rWy, + 2mV2(M + Ar')} + AK^r\M + Ar^Hrr - SAcoht^M + Ar')//^, = 

(t,r); 2r%M + Ar^H'^l + 2r\M + Ar^iM + 3Ar^)H^^ + r{M + Ar^) {-2M^ + Ar^x 
(2m2 - 5A) + Mr^{A + 2m^) + 2k^{M + Ar^) - 2r'^ujl] Ht, + Aik{M + Ar^ fHt^ 
+4AuJhr^Hu + 4AuJhr\M + Ar^fHrr = 



{t, 0); 2r2(M + Ar'^fHl - 2r{M + Ar'^fH',^ + + Ar''){2k'' + (A + 2m2)r2) 
-2r'a;^} + 4iA;r(M + Ar')'//*, + AAuhr^{M + Ar')i/,^ = 

(r,r); -2r^(M + Ar^)^//;^ - 2r=^(M + Ar2)3(A^ + TAr^)//;^ - r'(A^ + Ar^)^ {-AM^+ 
Mr^{2m^ + 5A) + Ar^(2m2 + 13A) + 2k'^{M + Ar^) - 2a;^r^} i/,,. + AA^r^Hu + 
-8Acu,,r^(A^ + Ar^)Htr + 4(A^ + Ar^H^^ - 8ikr{M + Ar^fHr^, = 



(r, 0); 2r^(>l + Ar^H';^ - 2r\M + Ar^iM - 3Ar')i/,^ + r{M + Ar^) {{M + Ar') 
x{2k^ - 6M) + r^(M + Ar^)(2m^ - 5A) - 20;^} Hr^, + Aikr^(M + Ar^fHrr 
+AAujhr^Ht^ + Aik(M + Ar'^fH^^ = 

0); 2r\M + Ar^fH'l^ - 2r{M + Ar2)(37W + Ar^)H'^^ + {{M + Ar^){2k^ + AM+ 
2mV2 + Ar^) - 2ujlr^} H^^ - Ar'^{M + Ar'^fH^ + 8i/cr(A^ + Ar'^fHr^, = 0. 
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Also the TT gauge condition of V^/i^i^ = and h = are given by 

t; r\M + Ar^K + <^ + Ar^){M + 3Ar^)Htr - ik{M + Ar^)Ht^ + r^UhHu = 
r; r^{M + Ar^fH'„ + r'^{M + Ar2)^(A^ + AAr'^)Hrr + Ar^if^ + (A^ + Ar'^fH^^ 

-ujhr^{M + Ar2)i7t, - ikr{M + Ar^Hr^ = 
0; -r2(7W + Ar2)2if;^ - r{M + Ar^){M + 3Ar^)Hr^ + ik{M + Ar^)/^^^ + uj^Ht^ = 

r'i/it - r^{M + Ar2)2if,, + (A^ + Ar')//^^ = 

For type I(odd) metric ansatz ( 173|) . i.e., iJtt = Htr = H„ = H^^ = 0, (t,r) equation becomes 

Mk{M + Ar^fHt^ = 0, 

and the solution is Ht^ = for k ^ 0. However, this corresponds to the null odd-solution 
because in this case, we obtain if^^ = from (t, 0) equation. Therefore, the s-mode {k = 0) 
solution is only admitted for type I metric ansatz (175]) . 

On the other hand, when focusing on {t, 0) equation and considering type II (even) metric 
ansatz with Ht^j) = Hr^ = 0, equation {t, 0) reduces to 

Aikr{M+ Ar'^fHtr = 0. 

Unless = 0, we obtain Htr = 0. In this case, Hu = is found from t component equation 
of V^h^y = 0. Furthermore, we have Hrr = and H^f,^ = from (t,t) and (r, r) equations. 
This corresponds to the null even-solution. So we also have to restrict to the s-mode case. In 
other words, the s-mode case leads to type I and II metric splitting for obtaining two modes 
of a massive graviton. 
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